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ABSTRACT 


The  effect  of  mobility  ratio  on  areal  sweep  efficiency  and  con- 
ductance ratio  has  been  made  for  direct  and  staggered  line  drive  well 
patterns.     This  study  was  made  using  miscible  fluids,    artificially 
consolidated  sandstone  models,    and  a  photographic  method  of  data 
recording.     The  mobility  ratio  was  varied  between  0.  1  and  12,    which 
is  the  range  common  to  water  flood  operations. 

A  method  of  predicting  the  performance  of  waterfloods  in 
stratified  reservoirs  using  conductance  ratios  is  treated.     This 
technique  is  not  limited  in  number  of  strata  and  indicates  one  method 
of  considering  a  free  gas  saturation. 
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INTRODUCTION 

The  importance  of  secondary  recovery  techniques  in  petroleum 
engineering  is  continually  increasing  and  is  taking  an  increasing  per- 
centage of  reservoir  research  activity.     Since  primary  reservoir  deple 
tion  normally  will  leave  approximately  75  per  cent  of  the  original  oil 
in  place,   the  reason  for  this  increasing  role  of  secondary  recovery  is 
obvious. 

Perhaps  the  oldest  and  most  used  method  of  recovering  oil  by 
secondary  techniques  is  the  waterflood.     This  involves  injection  of 
water  into  the  reservoir  to  displace  the  oil.     Water  injection  requires 
expenditure  of  energy  and  funds,    therefore,    reliable  methods  for  pre- 
dicting waterflood  results  are  necessary. 

Extensive  study  has  gone  into  conformance  factor  or  sweep 
efficiency  predictions  for  various  well  patterns  at  varying  mobility 
ratios.     This  permits  predictions  for  single  stratum  reservoirs. 

Most  reservoirs  are  composed  of  several  strata,   however,   and 
performance  predictions  for  stratified  reservoirs  are  not  feasible 
with  only  sweep  efficiency  data  available.     In  predictions  for  stratified 
systems,    it  is  necessary  to  place  performance  on  a  time  basis.     This 
requires  that  flow  potential  or  flow  rate  for  a  pressure  drop  be  avail- 
able for  the  system  under  investigation. 


Flow  potential  data  (conductance  ratios)  are  available  for  all 
regular  patterns  with  the  exception  of  the  line  drive.     The  pur- 

pose of  this  study  was  to  obtain  line  drive  conductance  ratios  and  to  test 
methods  for  approximating  conductance  ratios  analytically.     A  method 
of  predicting  performance  of  a  stratified  system  using  these  data  is 
also  presented. 


WATER  FLOOD  PREDICTIONS 

In  studies  of  petroleum  reservoirs,    much  qualitative  and  quanti- 
tative use  has  been  made  of  models  to  more  completely  understand  the 
phenomenon  of  oil  production.     Since  no  two  reservoirs  are  exactly  the 
same,    model  studies  are  limited  in  scope.     Therefore,    in  a  rigorous 
interpretation,    models  can  only  be  used  for  studies  of  a  particular 
field;  however,   by  generalizing  parameters,    study  of  a  reservoir  seg- 
ment will  yield  data  for  qualitative  application  in  a  number  of  situations. 

A  great  deal  of  work  has  gone  into  the  determination  of  scaling 

•*      •     4  a  i     «.  a-        3,4,10,17,23,29,31 

criteria  for  model  studies.  Probably  the  most 

complete  analysis  covering  general  situations  has  been  presented  by 

17 
Geertsma,    Croes,    and  Schwarz.  A  strict  interpretation  of  this 

analysis,   however,    leads  one  to  the  conclusion  that  the  model  must 

be  identical  to  the  prototype  in  every  respect.     Since  this  is  completely 

impractical,    it  is  necessary  to  decide  which  reservoir  factors  are  of 

greatest  importance  and  construct  the  model  accordingly.     A  detailed 

description  and  application  of  model  scaling  for  edgewater  drive  reser- 

4 
voirs  has  been  presented  by   Caudle. 

The  efficiency  of  a  secondary  recovery  project,    in  which  a  fluid 
(usually  water)  is  injected  to  displace  the  oil  in  the  reservoir,    can  be 
defined  by  the  volumetric  sweep  efficiency.     In  reservoir  studies, 
models  are  normally  two  dimensional  representations  and  the  displace- 
ment efficiency  may  be  defined  in  terms  of  areal  sweep  efficiency.     This 

3 


4 
is  the  ratio  of  the  area  swept  by  the  injected  fluid  to  the  total  model  area, 

Generally,   the  areal  sweep  efficiency  is  determined  by  the  well 
pattern  geometry  and  the  mobility  ratio.     This  assumes  that  the  reser- 
voir is  homogeneous  and  that  gravitational  forces  are  not  present.     The 
mobility  ratio  is  defined  as  the   ratio  of  the  sum  of  the  mobilities  of  the 
fluids  behind  the  front  to  the  sum  of  the  mobilities  of  the  fluids  ahead 
of  the  front.     The  mobility    (\)    of  a  fluid  in  a  porous  media  is  expressed 
as  the  effective  rock  permeability    (k)    to  the  fluid  of  interest  at  its 
saturation  divided  by  the  fluid  viscosity    (|i). 


(k/u)  swept  s  .»,.,.,«.  /,i 

M    =     »   ;    l K — r     =    : —    =    Mobility  Ratio  (1 

(k/p.)  unswept  \ 


Early  model  studies  utilized  the  electrolytic  model  to  represent 
waterfloods.     This  model  uses  the  principle  of  ionic  flow  to  simulate 

fluid  flow.     Use  of  this  model  is  limited  to  mobility  ratios  of  one. 

27 
Muskat,    Wyckoff,   and  Botset       were  among  the  first  to  utilize  this 

method  to  determine  breakthrough  sweep  efficiencies  for  various  well 

26 
patterns.     In  a  later  publication,    Muskat  and  Wyckoff       analyzed  the 

27 
previous  electrolytic  model  results        and  concluded  that  the  shape  of 

the  flood  front  and  the  breakthrough  areal  sweep  efficiency  were  depend- 
ent on  the  geometry  of  the  system  and  not  on  pressure  differential,   well 

24, 25, 26 
spacing,    or  rock  properties.     It  is  noted  that  Muskat  indicates 


that  the  results  obtained  in  the  electrolytic  models  do  not  consider  the 
effects  of  two  fluids  with  different  viscosities.     In  other  words,    the  elec 
trolytic  model  is  valid  for  a  mobility  ratio  of  one  only. 

The  electrolytic  model  and  the  similar  but  more  sophisticated 
potentiometric  model  have  since  been  used  in  the  study  of  various  reser 

voir  conditions.     Limited  use  has  also  been  made  of  heat  flow  models  to 

22 
simulate  fluid  flow  with  a  mobility  ratio  of  one. 

Sweep  Efficiency 

Early  studies  on  the  effect  of  mobility  ratio  on  areal  sweep  effi- 
ciency utilized  potentiometric  models.     Aronofsky     reported  results  of 
a  potentiometric  model  study  of  mobility  ratio  effects.     A  stepwise  pro- 
cedure was  utilized  to  vary  the  ratio  of  the  depth  of  the  electrolyte 
behind  the  front  to  that  ahead  of  the  front.      This  variation  in  depth  made 
the  electrical  conductivity  ahead  of  the  front  different  from  that  behind 
the  front,    thereby  simulating  fluids  of  different  mobilities  flowing  in  a 
porous  media.     After  computing  a  new  frontal  position,    the  electrolyte 
depth  was  changed  to  fit  the  new  frontal  position  and  the  process  was 
repeated.     In  this  iterative  manner,    sweep  efficiencies  for  mobility 
ratios  other  than  one  were  approximated  by  advancing  the  front  in  a 
series  of  steady-state  steps. 

The  effect  of  mobility  ratio  on  areal  sweep  efficiency  has  also 

28 
been  studied  with  an  electric  analog  model  by  Nobles  and  Janzen. 


This  study  utilized  a  resistance  network  and  an  iterative  process  similar 

to  Aronofsky's     to  simulate  frontal  movement. 

8 
Cheek  and  Menzie     studied  the  effects  of  mobility  ratio  on  areal 

sweep  through  the  use  of  a  fluid  mapper.     This  model  uses  fluid  flow 

between  two  closely  spaced  parallel  plates  to  simulate  fluid  flow  in  a 

porous  media.     Results  are  reported  for  five  spot  and  direct  line  drive 

well  patterns. 

32 
Slobod  and  Caudle        introduced  the  X-ray  shadowgraph  technique 

to  study  fluid  flow  phenomenon.     This  new  technique  used  a  porous  plate 
and  fluid  flow  with  one  of  the  fluids  containing  an  X-ray  absorbing  mate- 
rial so  that  a  contrast  between  displaced  and  displacing  fluid  could  be 
observed  on  X-ray  shadowgraphs  taken  during  the  experiment.     By 
taking  successive  X-ray  exposures,    a  complete  history  of  the  flood  was 

obtained. 

13 
Dyes,    Caudle,   and  Erickson       used  the  X-ray  shadowgraph 

method  to  illustrate  that  a  significant  amount  of  production  could  take 

place  after  initial  breakthrough.     The  results  of  this  study  demonstrated 

that,   with  most  water  floods,    an  ultimate  areal  sweep  efficiency  of  85- 

100  per  cent  could  be  obtained  for  various  mobility  ratios  at  an  economic 

water-oil  ratio.      The  only  previous  reports  which  indicated  consideration 

21 
of  production  after  breakthrough  were  made  by  Hurst       for  a  mobility 

15 
ratio  of  one  and  Fay  and  Prats        for  a  mobility  ratio  of  4.  0  using 

numerical  techniques. 


The  X-ray  shadowgraph  has  been  used  extensively  to  study  a 

5,6,7,11,13,14,16 
variety  of  reservoir  conditions  at  varying  mobility  ratios. 

Included  among  these  studies  are  five -spot  pattern  floods,    stratified 

five -spot  floods,    and  the  effect  of  fractures  on  flood  performance. 

1 8 
Habermann        utilized  a  method  similar  to  the  X-ray  shadow- 
graph for  studying  fluid  flow.     This  method  consisted  of  artificially 
consolidated  sand  models  covered  with  clear  plastic.     One  of  the  fluids, 
either  displaced  or  displacing,    contained  a  dye  which  permitted  visual 
observation  of  the  front  as  it  progressed.     Photographs  taken  during 
the  flood  permitted  calculation  of  areal  sweep  efficiency. 

Conductance  Ratio 

The  use  of  areal  sweep  efficiency  data  permits  ultimate  re- 
covery predictions  for  single  stratum  reservoirs.     In  order  to  predict 
performance  of  stratified  reservoirs  and  also  to  predict  recovery  of 
single  stratum  reservoirs  on  a  time  basis,    it  is  necessary  to  know  the 
change  in  the  flow  potential  or  flow  rate  per  unit  pressure  drop  for  the 

system  under  study. 

24 
Prediction  of  flow  potential  was  first  discussed  by  Muskat        in 

terms  of  conductivity  (production  rate  per  well  per  unit-pressure 

differential).     This  permitted  predictions  for  systems  with  mobility 

ratios  of  one. 

2 
Aronofsky  and  Ramey     first  introduced  the  concept  of  conductance 


8 
ratio  in  a  potentiometric  study  of  mobility  ratio  effects  on  a  five -spot 
water  flood.     The  conductance  ratio    (y)    is  defined  as  the  flow  rate     (q; 
per  unit  pressure  drop    (AP)    between  the  injection  and  production  wells 
at  a  specified  time  to  that  conductivity  with  only  the  original  fluid  pre- 
sent. 


(q/AP)  time  t  _       ,  _  ,_% 

=    7*  .    _.   . — — — r    -    Conductance  Ratio  (2) 

(q/AP)  initial  v    ' 


Analytical  methods  of  studying  the  change  in  flow  potential  has 

been  presented  in  terms  of  dimensionless  injectivity  by  various  methods. 

30 
Prats,   Matthews,   Jewett,    and  Baker       have  computed  injectivity  for 

five -spot  systems  by  considering  that  the  reservoir  consists  of  several 

19 
layers,    each  having  a  uniform  but  different  permeability.     Hauber 

presented  another  approximate  method  which  considers  that  flow  uti- 
lizes the  stream  tubes  and  iso-potential  lines  obtained  by  the  solution 
of  differential  equations  for  a  mobility  ratio  of  one. 

Deppe    '"  presented  a  method  to  approximate  conductivity  for  vary- 
ing mobility  ratios  and  for  varying  well  patterns.     In  the  five -spot 

12 
pattern  used  as  an  example,    Deppe    '"  essentially  divided  the  pattern 

into  two  regions  of  radial  flow.      The  first  region  was  around  the  injec- 
tion well  and  covered  almost  one -half  of  the  area.     The  second  radial 
flow  region  was  the  same  except  that  it  was  around  the  production  well. 


12 
The  results  of  Deppe's        study  provided  injectivity  (conductance 

ratio)  curves  which  were  not  too  unreasonable  for  five -spot  and  nine- 
spot  well  patterns.     Closer  agreement  was  obtained  for  unfavorable 
mobility  ratios,   which  is  normally  true  when  using  this  type  of  approxi- 
mation methods. 

In  an  attempt  to  improve  on  the  work  of  Deppe,    this  study  has 
divided  the  pattern  into  3  regions  as  illustrated  by  Figure  1  for  a 
staggered  line  drive  pattern.     The  first  region  surrounds  the  injection 
well  where  the  largest  pressure  drop  is  encountered.     This  region  is 
considered  to  have  radial  flow  and  has  a  radius  equal  to  one -fourth  the 
well  row  spacing.     The  second  region  is  considered  to  be  linear,   with 
only  slight  pressure  change,   and  extends  to  the  areal  sweep  at  break- 
through.    The  third  region  is  one  of  radial  flow  around  the  production 
well  where  the  pressure  drop  is  considered  to  be  a  function  of  frac- 
tional flow  from  the  swept  region. 

Conductance  ratio  for  the  system  is  computed  by  considering 
a  constant  flow  rate  which  reduces  the  equation  to 

V    =    *Ii  (3) 

ftPt 

Pressure  drop  across  the  system  is  obtained  by  summing  the  pressure 
drops  in  the  three  regions  at  the  time  of  interest. 


10 
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&  Angle  open  to  flow  from  swept  region  after  breakthrough 


FIGURE  1.     STAGGERED  LINE  DRIVE  --  ASSUMED  FLOW  REGIONS 


11 

With  only  single  fluid  flow  in  the  system  (or  at  initial  conditions) 
the  pressure  drop  in  each  region  can  be  represented  by  modifications  of 
Darcy's  Law  as 


q  In  ri/rw  . 

APAi \    7.08h  (4) 

u 


-    q  l/w  ... 

APBi    ~    \  1.127h  (5) 

u 


_     q  In  rp/rw  .,. 

APCi    "       \    ?.08h  (b) 

u 


which  when  added,    gives  an  initial  pressure  drop  for  the  system  of 

APi    =    \   7%8h     ^ln  ri/rw  +  6*28  L/w  +  ln  rP/rw)  (7) 

u 

After  the  flood  has  started  and  two  fluids  are  flowing  in  the  model,   the 
pressure  drop  can  be  represented  as  a  summation  of  the  pressure  drop 
in  each  fluid  region.     When  the  front  is  in  the  first  region  the  pressure 
drop  therein  can  be  represented  as 


q  ln  rf  /  rw  q  ln  ri/rf 

APAt  \    7.08h       +      X   7.08h 


s  u 

(8) 

q  (1/M)  ln  rf/rw  +  ln  ri/rf 

X    7.08h 
u 


12 
When  the  front  is  in  the  second  region  the  pressure  drop  therein  is 


a-p  qi  ,       q(  L-  A) 

A*Bt    "    \_1.  127hw  \.l.  127h 

(9) 


rTEw 

s  u 


Sr-] 


i/w  +  L/w 


"    q  X    1.127h 

u 


After  breakthrough,   the  flow  rate  in  the  third  region  must  be  summed 
for  parallel  flow  since  two  fluids  are  flowing  into  the  producing  well. 
Consider  that  portion  of  the  well  illustrated  in  Figure  1  or  one-quarter 
of  the  well.     Let    a    equal  the  angle  open  to  flow  from  the  swept  region 

and    (rr/2-a)    equal  the  angle  open  to  flow  from  the  unswept  region. 

30 
Then  the  fractional  flow  from  the  swept  region  is 


\  a 
fs    =    _ f    ,    p, ■  (10) 


then 


fs  rr/2  .,  , . 

a    =    M-fs(M-l)  (11) 


13 


The  pressure  drop  in  the  swept  and  unswept  region  is 


fs  q  In  rp/rw  rr/2 
a 


APCs 7.08h\    - (12) 


(1-fs)  q  In  rp/rw  tt/2  ,.  ~. 

APCu ?.08h\     (n/2-a) (13) 

u 


For  parallel  flow,   the  pressure  drop  given  by  equation  (12)  must  equal 
that  of  equation  (13).     Therefore,   the  pressure  drop  in  the  third  region 
may  be  obtained  from  either  the  equation  for  the  swept  or  the  unswept 
portion.     By  substituting  equation  (11)  into  (13)  we  obtain 

(1-fs)  q  In  rp/rw 


APCt    "    7.08h  Xu  ri-fs/[M-fs(M-l)]j]  (14) 

In  summing  equations  (8),    (9),    and  (14)  the  following  conditions 
apply: 

1.  When  Es<Es-     jj/w  =  0,    fs  =  0,    and  rf  =  f(Es). 

2.  When  Es     <Es<Es        :     rf  =  ri,    fs  =  0,    and  jfc/w  =  f(Es). 

3.  When  Es>Esnrp:     rf  =  ri,    £/w  =  L/w,    and  fs  =  f  (Es). 

ri  1 


The  total  pressure  drop  at  any  time  is  then 
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AP      =    3 

a±t  X   7.08h 


u 


+    6.  28 


w-1 


(1/M)  In  rf/rw  +  In  ri/rf 


4/w  +  6.  28  L/w 


(1-fs)  In  rp/rw 
+    l-*s/|M-fs(M-l)] 


(15) 


The  conductance  ratio  from  equation  (3)  becomes  the  ratio  of  equation 
(7)  to  equation  (15)  or 


Y    =    APi/APt    =    (In  ri/rw  +  6.  28  L/w  +  In  rp/rw)/ 


(1/M)  In  rf/rw  +  In  ri/rf  +  6.  28 


a-0 


4/w 


♦  6-«>Ww  +  r^TT^fe) 


(16) 


Values  for     L/w    were  obtained  by  using  both  arbitrary  values  and  by 
assuming  that  the  pressure  drop  across  a  linear  system  is  equal  to 
that  across  a  line  drive  with  well  radius  equal  to    ri.     Muskat's  equation 
for  flow  rate  in  a  line  drive  is: 


AP    = 


qii    d/q  -   1.  17  +  0.637  log  q/ri 


kH 


2.254 


(17) 
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Setting  equation  (17)  equal  to  (5)  and  solving  for  L/w  yields 

L/w    =    d/a-1.17+026371oga/ri  (lg) 

which  is  valid  when    d/a>l.      Calculations  were  performed  by  programming 
the  above  for  the  Controlled  Data  Corporation  1604  digital  computer. 
Appendix  D  contains  the  flow  diagram  and  FORTRAN  language  program 
utilized  in  these  computations.     A  conversion  table  for  the  FORTRAN 
symbols  is  also  provided  in  Appendix  D. 

In  the  calculations  for  conductance  ratios,   various  values  of  L/w 
were  assumed  in  an  attempt  to  match  the  model  study  results.     Using 
equation  (18),    L/w  was  0.357.     Other  values  used  for  L/w  were  0  and 
0.  1.     As  can  be  seen  from  the  plots  of  Appendix  D,   varying  L/w  made 
only  small  changes  in  the  conductance  ratio.     Comparison  plots  for  L/w 

of  0.  1  and  as  obtained  from  the  model  studies,    indicate  results  similar 

12 
to  those  obtained  by  Deppe.  Agreement  is  fair  for  unfavorable  mobility 

ratios  and  not  too  good  for  favorable  mobility  ratios. 

This  method  can  be  used  for  rough  approximation  of  conductance 
ratios  for  varying  well  patterns  by  modifying  the  initial  values  in  the 
program  to  correspond  to  those  of  the  pattern  under  study.     However, 
it  is  believed  that  fluid  flow  model  study  results  would  be  more  reliable. 

The  effect  of  mobility  ratio  on  conductance  ratio  has  been 
studied  for  various  well  patterns  using  fluid  flow  models.      Caudle  and 
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7 
Witte     have  published  data  for  the  five -spot,   the  inverted  skewed  seven- 

20 
spot  has  been  studied  by  Hickman,        and  the  inverted  nine -spot  was 

34 
studied  by  Watson.  This  study  provides  data  for  line  drive  well 

patterns. 

Stratified  Reservoirs 

In  a  reservoir  composed  of  strata  of  varying  permeabilities 
and  porosities,   performance  prediction  becomes  more  complex  than 
for  a  single  stratum.     To  illustrate,    consider  a  two  strata  system  of 
equal  characteristics  with  fluids  of  equal  viscosity  and  with  the  perm- 
eability in  the  first  stratum  twice  that  of  the  second.     From  Darcy's 
Law 

it  can  readily  be  seen  that  the  first  stratum  will  have  a  frontal  advance 
twice  that  of  the  second  strata  and  that  at  breakthrough  for  the  first 
stratum  the  second  stratum  will  have  a  volumetric  sweep  of  only  half 
of  the  breakthrough  sweep  efficiency.     In  this  simple  example,   a  mis- 
cible  system  was  assumed  in  order  that  mobility  ratios  of  one  could 
be  used  for  each  strata.     For  mobility  ratios  other  than  one,    the  linear 
relationship  of  flow  velocity  to  permeability  no  longer  exists  in  a 
stratified  system  and  then  becomes  a  function  of  distance  traveled  or 
area  swept. 
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In  order  to  apply  the  data  obtained  in  this  study  to  a  stratified 
reservoir,    it  is  necessary  to  put  the  performance  of  each  strata  on  an 
equal  time  basis.     This  is  most  easily  done  through  the  use  of  the  con- 
ductance ratios  obtained  in  this  and  other  studies.     It  is  noted  that  this 
method  is  developed  and  accurate  for  separated  strata  or  for  strata 
without  cross  flow.     When  cross  flow  between  strata  exists  and  the 
mobility  ratio  is  greater  than  one,    this  method  provides  a  good 
approximation  for  predicting  water  flood  results.     A  similar  method 

using  theoretical  dimensionless  injectivities  has  been  presented  by 

30 
Pratts,    et  al.  The  method  developed  herein  permits  performance 

prediction  as  accurately  but  more  simply  than  the  method  described 

30 
by  Prats,   et  al. 

To  develop  this  time  relationship,    consider  Darcy's  Law  (19) 

and  let 

v     =    q/A 

and 

A    =    Vd/A 

Darcy's  Law  can  then  be  written  in  terms  of  conductance  as 

q/AP    = 7--TT    =    conductance  (20) 

u  i  hi 

Consider  now  two  adjacent  parallel  strata  in  a  reservoir  with  single 
fluid  flowing;  the  variables  between  these  two  strata  are  permeability  (k) 
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and  displaceable  volume     (Vd). 

Therefore  the  conductance  between  a  stratum  1  and  any  stratum 
j    will  vary  as  their  permeability  and  displaceable  volume  ratios  vary 
or 


i.  k.Vd.         q. 


J 


(21) 


By  definition, 


Y    =     (q/AP)t  /   (q/AP). 


(22) 


then 


k.Vd. 


h/ipj]t  ■  yj  «X 


AP:J 


(23) 


and 


dWi 
qj    =    "dT 


k.Vd. 


LiPjJi 


(24) 


solving  for    dt,     integrating,   and  letting    c    equal    AP= 


r^P 


.APlJ 


ct     = 


klvdl 

k.Vd. 

J     J 


1/v.    dWi 


(25) 
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Other  equations  needed  to  express  stratified  reservoir  perform- 
ance are  for  total  flow  from  the  swept  region    (fs)    and  total  conductance 
ratio    (y). 

By  definition, 

fg     =    qs/q    =    E    qs.     /     Y    q.  (26) 

then 

qs.     =     fs.     q.  (27) 


and 


APi 


from  (7) 


qj  ■  yj  Cq3]i  rAPjT;  (28) 


k.Vd. 

CqA   =  At^    Cq,i4  (29) 


lj  'i  -  tt-w;    L4i 


i  T~i 


then  by  substituting  into  (26)  we  have 

£  y-  fs.  k.  Vd. 

fs  =  .A'vl  '  <30> 

The  expression  for  overall  conductance  ratio  can  be  derived  in  this 

same  manner  as 

T  Y-  k.  Vd. 

v  -  5  L  k  J  <31> 

J       J 
An  example  performance  prediction  for  a  stratified  reservoir 

using  this  method  is  described  in  Appendix  E.     A  possible  method  of 

considering  an  initial  free  gas  saturation  is  also  mentioned  in  Appendix 

E. 


EXPERIMENTAL  APPARATUS 

Model  Theory 

The  majority  of  oil  reservoirs  are  drilled  with  a  regular  pattern 
of  wells.     Probably  the  most  common  pattern  utilized  in  this  country  is 
a  square  array  of  wells  on  a  40  acre  spacing.     This  well  pattern,    illus- 
trated by  Figure  2,    readily  lends  itself  to  secondary  recovery  operations 
normally  referred  to  as  the  Direct  Line  Drive,   which  is  obtained  by  con- 
verting alternate  rows  of  wells  to  injection  wells. 

The  Staggered  Line  Drive  is  also  obtained  by  converting  alternate 
rows  of  wells  to  injection  wells;  however,    the  wells  are  not  drilled  in  a 
square  array  but  rather  in  an  offset  fashion.     The  well  pattern  for  the 
Staggered  Line  Drive  system  is  shown  by  Figure  3. 

With  the  regular  well  systems  shown  by  Figures  2  and  3,   it  can 
readily  be  seen  that  there  will  exist  numerous  mirrow  images  divided 
by  lines  of  symmetry.     These  lines  or  planes  of  symmetry  represent 

invarient  streamlines  or  boundaries  across  which  no  flow  will  occur, 

9 
as  long  as  production  and  injection  quantities  are  uniform.        Since  the 

planes  of  symmetry  are  in  effect  flow  boundaries  in  a  homogeneous 
system,    it  is  apparent  that  only  the  smallest  segment  bounded  by 
symmetry  planes  need  be  considered.     Figures  4  and  5  show  enlarged 
views  of  the  model  and  its   relationship  to  the  line  drive  systems  studied. 
As  can  be  seen,   the  model  utilized  represents  a  line  drive  pattern  where- 
in the  row  spacing  equals  the  well  spacing. 
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FIGURE  2.     DIRECT  LINE  DRIVE  WELL  ARRAY 
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FIGURE  3.     STAGGERED  LINE  DRIVE  WELL  ARRAY 
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24 


o  — ■ 

* 


-  a 


o 


o —   —    —    0 


o  - 


-  o  — 


O        Production  Well 

0        Injection  Well 
— Lines  of  Symmetry 

Model  Used 
d  Distance  Between  Rows 

a  Distance  Between  Wells 


FIGURE  5.     STAGGERED  LINE  DRIVE  PATTERN 
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Model  Construction 

The  model  utilized  in  this  study  was  an  artificially  consolidated 
sand  plate  constructed  generally  as  described  by  Watson       and  Hickman. 

The  artificial  sandstone  was  composed  of  a  clear  epoxy  resin  and 

60-100  mesh  Ottawa  sand.     Actual  construction  of  the  artificial  sandstone 

34  20  4 

has  been  reported  in  detail  by  Watson,        Hickman,        and  Caudle.        In 

order  to  avoid  any  possible  loose  streaks  which  may  have  occurred  during 

packing  of  the  artificial  sand  a  minimum  of  2  inches  was  cut  from  the 

model  sides  by  using  a  diamond  saw.     To  seal  the  model,    an  extremely 

viscous  epoxy  resin  was  spread  over  all  exposed  sand  faces.     The  high 

viscosity  of  the  epoxy  resin  was  necessary  to  prevent  imbibition  of  the 

coating  into  the  model  and  was  obtained  by  applying  the  coating  just 

prior  to  it  being  completely  "set-up.  " 

Injection  and  producing  "wells"  were  constructed  from  1/2  inch 
thick  acrylic  plastic  plate  as  shown  by  Figure  6.     Swagelok  fittings  for 
1/8  inch  tubing  were  utilized  to  connect  the  "wells"  to  the  injection 
systems.     Prior  to  attaching  the  plastic  plates  to  the  model,   well  bores 
were  reamed  out  at  the  well  corners. 

As  seen  from  Figure  6,  the  model  utilized  was  32  cm  in  length, 
16  cm  in  width,  0.  6  cm  thick,  and  had  a  well  bore  diameter  of  0.  2  cm. 
This    represents  a  d/rw  ratio  of  320. 
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1/8"  Fitting 


0.6  cm  Consolidated 
Sandstone 


FIGURE  6.     LINE  DRIVE  MODEL 
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Fluids 

The  fluids  utilized  in  this  study  were  miscible  in  all  proportions. 
The  use  of  miscible  fluids  to  simulate  reservoir  fluids  in  model  studies 
of  mobility  ratio  effects  has  been  well  established  with  experimental  data 
by  Slobod  and  Caudle    '"  and  in  other  publications.     '     '     '       '       ' 

By  using  miscible  fluids  in  this  study,    there  are  no  effects  of 

capillary  forces  as  will  be  found  in  immiscible  systems.     The  absence 

of  capillary  effects  eliminates  the  problem  of  saturation  gradients  in 

determining  the  effective  mobility  ratio.     Also  by  utilizing  miscible 

fluids,   the  effective  permeability  to  each  fluid  is  the  same,   thereby 

eliminating  the  need  to  determine  relative  permeabilities.     This  system 

therefore  models  a  water  flood  with  little  or  no  flow  of  oil  behind  the 

water  front.     The  mobility  ratio,      ,,  ,''. *- ,   then  reduces  to  the 

(k/^i)  unswept 

viscosity  ratio  of    4*-r -*—  .     Therefore,   only  the  viscosity  of  the 

'  {\l)  swept  ' 

fluids  needed  to  be  changed  in  order  to  vary  the  mobility  ratio  of  the 
system. 

To  obtain  fluids  of  the  desired  viscosities,   a  mixture  of  gly- 
cerine and  water  was  used.     The  quantity  of  glycerine  added  for  the 
desired  viscosity  was  determined  on  a  weight  per  cent  basis  through 
the  use  of  standard  chemical  tables  and  checked  with  an  Ostwald 
viscosimeter. 

In  the  photographic  technique  used  in  this  study,    one  of  the  fluids 
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contained  a  dye  to  establish  a  contrast  between  the  displacing  and  the 
displaced  fluids.     Since  the  epoxy  used  to  consolidate  the  sand  renders 
the  grains  oil  wet,    oil  soluble  dyes  could  not  be  used.     It  was  found  that 
an  ordinary  household  food  coloring  used  in  about  5  per  cent  volumetric 
proportions  with  the  water  obtained  the  desired  contrast. 

Injection  System 

20 
The  constant  rate  pumps,    utilized  and  described  by  Hickman 

34 
and  Watson       were  used  in  this  study.     This  pump  is  basically  a  hydraulic 

cylinder  using  a  screw  type  mechanism  to  move  the  piston.     It  is  driven 
by  a  gear  and  chain  arrangement  from  an  electric  motor  through  a 
variable  speed  reduction  gear.     By  changing  the  gears  in  the  drive  mech- 
anism and/or  the  speed  output  of  the  reduction  gear,   the  desired  injection 
rate  could  be  obtained. 

A  pressure  gauge  was  located  at  the  injection  well.     The  producing 
well  was  permitted  to  discharge  at  atmospheric  pressure,   thereby  per- 
mitting the  gauge  to  indicate  the  pressure  drop  across  the  model. 

Data  Recording 

A  35  mm  camera  was  used  to  record  the  frontal  position,   pressure 
reading,  and  time  during  the  runs.     The  camera  was  actuated  by  a  D.  C. 
timing  motor  which  advanced  the  film  and  made  one  exposure  per  minute. 
Exposures  could  be  made  at  longer  intervals  by  using  a  "sign  flasher"  to 
interrupt  current  flow  to  the  timing  motor.     This  was  done  during  the 
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slower  runs.     A  red  filter  was  used  to  obtain  better  contrast  between 
the  blue  and  clear  fluids  on  black  and  white  medium  speed  panchro- 
matic film.     A  more  complete  description  of  the  photographic  and 

4 
developing  technique  has  been  previously  written  by  Caudle. 

Figure  7  is  a  schematic  diagram  of  the  equipment  used  in  this 
study. 

Typical  photographs  are  shown  in  Figures  8  and  9.     Figure  8 
is  an  exposure  of  a  direct  line  drive  at  a  mobility  ratio  of  one  and  Figure 
9  is  for  a  staggered  line  drive  at  a  mobility  ratio  of  twelve.     Both  expo- 
sures indicate  frontal  positions  shortly  before  breakthrough. 
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FIGURE  8.     PHOTOGRAPH  OF  MODEL  DURING  A  DIRECT  LINE 
DRIVE  RUN,     M  =  1. 
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FIGURE  9.     PHOTOGRAPH  OF  MODEL  DURING  A  STAGGERED 
LINE  DRIVE  RUN,     M  =  12. 


EXPERIMENTAL  PROCEDURE 

The  performance  of  both  direct  and  staggered  line  drive  patterns 
were  studied  for  mobility  ratios  between  0.  land  12  to  determine  the 
effect  on  areal  sweep  efficiency,   fractional  flow,   and  conductance  ratio. 

The  model  was  initially  evacuated  by  means  of  a  vacuum  pump. 
After  complete  evacuation,   the  model  was  saturated  with  water.     At 
this  time,   all  lines  and  the  pressure  gauge  were  connected  to  the  model 
and  the  injection  pump  (Figure  7).     The  injection  pump  was  loaded  with 
dyed  water  as  the  initial  displacing  fluid,   the  model,    pressure  gauge, 
and  clock  were  placed  on  a  light  table  beneath  the  35  mm  camera  and 
the  direct  line  drive  runs  were  started. 

As  injected  fluid  moved  into  the  model,   the  automatic  camera 
recorded  the  advancing  front  positions  and  injection  pressures  at  regular 
intervals.     The  camera  was  stopped  after  approximately  95  per  cent  of 
the  model  had  been  swept.     In  either  case,   injection  was  continued  until 
complete  displacement  of  the  "oil  analog"  had  occurred. 

In  the  second  run,   the  injection  system  was  loaded  with  fluid  of 
higher  viscosity  and  the  above  process  was  repeated  with  the  clear 
viscous  fluid  displacing  the  dyed  fluid.     By  alternating  injection  between 
clear  fluid  of  varying  viscosity  and  the  dyed  fluid,   the  range  of  mobility 
ratios  from  0.  1  to  12.  0  was  covered  in  five  runs. 

After  completion  of  a  set  of  runs  covering  the  range  of  mobility 
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ratios  under  study,    the  roles  of  the  injection  and  producing  wells  were 
reversed  and  the  five  runs  were  repeated.      This  procedure  was  then 
followed  for  the  staggered  line  drive  system.     The  two  runs  at  the 
same  conditions  were  made  in  opposite  directions  to  average  out  any 
difference  in  the  pressure  drop  around  the  individual  wells. 

The  flow  rate  was  set  prior  to  the  start  of  each  run  and  was 
maintained  constant  throughout  the  run.     The  rate  was  chosen  so  that 
the  injection  pressure  would  not  exceed  the  gauge  scale  and  so  that  a 
sufficient  number  of  photographic  exposures  would  be  made  during  the 
run  to  permit  interpretation.     The  areal  sweep  efficiency  is  independent 
of  flow  rate  in  the  line  drive  floods.     Since  the  conductance  ratio 
(Equation  2)  is  a  dimensionless  ratio,    the  actual  rate  utilized  is  of  no 
concern.     It  is  important,   however,    to  ensure  that  the  flow  rate  remains 
constant  during  each  run  so  that  the  time  scale  can  later  be  converted  to 
pore  volumes  of  water  injected  and  the  conductance  ratio  computed  on 
the  ratio  of  pressure  drop. 


ANALYSIS  OF  DATA 

Areal  Sweep  Efficiency 

All  data  was  recorded  on  film  during  each  run.     After  developing 
the  film,   a  microfilm  reader  was  used  to  enlarge  the  recorded  image. 
Tracings  were  then  made  of  the  model  and  the  frontal  positions  during 
sweepout.     Appendix  A  contains  typical  tracings  of  frontal  positions  at 
various  mobility  ratios  for  both  the  direct  line  drive  and  the  staggered 
line  drive.     The  area  of  the  model  and  the  area  behind  the  frontal  posi- 
tion were  obtained  with  a  polar  planimeter.     The  areal  sweep  efficiency 
was  then  obtained  for  each  frontal  position  by  dividing  the  area  swept  by 
the  area  of  the  model.     Tables  of  data  are  contained  in  Appendix  B. 

Areal  sweep  efficiencies  were  plotted  versus  an  elapsed  time 
scale.     Since  runs  at  the  same  mobility  ratio,   but  in  opposite  directions, 
were  not  necessarily  made  at  the  same  rate,   a  time  correction  factor 
was  obtained  by  plotting  areal  sweep  data  for  the  second  run  on  the 
first  curve  obtained  and  then  noting  the  time  correction  necessary  which 
would  make  the  curves  most  closely  coincide.     This  permitted  easy 
correlation  of  data  for  runs  of  the  same  mobility  ratio  made  in  opposite 
directions. 

Since  each  run  was  made  at  a  constant  rate,  it  is  easily  seen  that 
the  time  scale  of  the  basic  data  plots  is  directly  related  to  the  volume  of 
displacing  fluid  injected  by  a  constant  ratio  for  each  particular  rate. 
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Prior  to  breakthrough  of  injected  fluid,    a  plot  of  areal  sweep  versus 
volume  injected  is  a  straight  line.     By  noting  the  area  swept  and  the 
elapsed  time  at  breakthrough,   all  subsequent  time  values  may  be  con- 
verted to  displaceable  volumes  of  injected  fluid  by  the  ratio  of  area 
swept  to  elapsed  time  at  breakthrough.     Appendix  B  contains  plots  of 
all  mobility  ratios  studied  for  the  two  well  patterns  considered,   with 
tables  of  data  derived  from  these  plots. 

Plots  of  areal  sweep  versus  displaceable  volumes  injected  can 
be  found  in  Appendix  C.     The  displaceable  volume  injected  is  a  more 

useful  term  than  actual  volume  (or  time  at  constant  injection  rate)  due 

13 

to  its  dimensionless  nature.     Dyes,    et  al       has  defined  displaceable 

volume  as  the  product  of  the  hydrocarbon  pore  volume  of  the  unit 
pattern  and  the  volumetric  displacement  efficiency. 

It  is  noted  that  areal  sweep  data  from  the  first  5  runs  did  not 
correspond  to  other  runs  for  the  direct  line  drive  made  in  the  opposite 
direction  as  regards  areal  sweep  efficiency.     After  study  of  the  photo- 
graphed flood  front  shapes,   it  was  noted  that  an  air  pocket  had  been 
present  near  the  production  well  during  these  runs,    thereby  distorting 
streamlines  and  causing  invalid  areal  sweep  data.     This  error  was 
noted  after  completion  of  all  runs  and  an  attempt  was  then  made  to  re- 
make these  runs.     However,    the  injection  well  became  "plugged"  during 
these  "reruns"  and  areal  sweep  data  for  the  direct  line  drive  for  mobility 
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ratios  of  3.  13,    0.  106,   and  12.  1  is  based  on  one  run  only.     The  cause 
of  this  "well  plugging"  is  attributed  to  impurities  in  the  injected  fluid, 
which  could  have  been  prevented  by  installing  a  filter  on  the  injection 
line.     This  air  pocket  did  not  interfere  with  runs  subsequent  to  run  5 
since  a  number  of  days  elapsed  between  run  5  and  6  and  the  model  was 
routinely  evacuated  before  making  run  6. 

Conductance  Ratio 

2 

The  conductance  ratio  is  defined  by  Aronof  sky  and  Ramey 

(q/AP).    /  (q/AP).    which  reduces  to    Ap.  /  AP      when  the  flow  rate  is 
maintained  constant.     Pressure  data  taken  from  the  injection  well  pro- 
vides   AP    at  a  particular  areal  sweep  efficiency  and  permits  direct 
calculation  of  the  conductance  ratio. 

From  the  above  equation,    it  can  be  seen  that  if  the  pressure 
increases  during  a  run,   the  conductance  ratio  will  decrease  and  a  de- 
creasing pressure  will  cause  an  increasing  conductance  ratio.     The 
conductance  ratio  decreases  if  the  mobility  ratio  is  favorable  and  in- 
creases if  the  mobility  ratio  is  unfavorable.     At  a  mobility  ratio  of  one, 
the  conductance  ratio  will  remain  constant  during  a  flood. 

Plots  of  conductance  ratio  versus  areal  sweep  can  be  found  in 
Appendix  C. 

The  ratio  of  the  distance  between  injection  and  production  wells 

2 

to  the  wellbore  radius    d/rw    has  been  shown  by  Aronof  sky  and  Ramey 
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to  effect  the  conductance  ratio.     Although  the  model  was  constructed  with 
a    d/rw    ratio  of  320,    slight  differences  in  the  well  bore  radius  and  plugging 
around  the  well  bore  will  tend  to  change  this  ratio. 

The  steady  state  flow  equation  for  single  phase  flow  in  line  drive 

24  12 

patterns  was  derived  by  Muskat       and  presented  by  Deppe        in  practical 

units  as 

_   0.001538  k/tih  AP 

q     '    0.682  d/a    -0.798  +  In  d/rw  (      } 

This  equation  is  valid  for  both  direct  and  staggered  line  drive  well 
patterns  with    d/a>l. 

Using  the  information  obtained  from  the  model  runs,    the  effec- 
tive   d/rw    ratio  for  the  direct  line  drive  was  found  to  be  922.     The 
effective    d/rw    ratio  for  the  staggered  line  drive  was  calculated  as   1110. 
These  values  for    d/rw    are  somewhat  smaller  than  those  normally  en- 
countered in  the  field.     However,    it  is  considered  that  reasonable 
accuracy  will  be  obtained  by  using  the  conductance  ratios  obtained  in  this 
study  without  correction. 

If  it  is  desired  to  correct  the  conductance  ratios  obtained  in  this 
study  to  the  correct  value  for  other    d/rw    ratios,    using  the  method  pre- 
sented  by  Aronofsky  and  Ramey     for  line  drives  yields  the  following  re- 
lationship: 
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(rw) 
0.68Z(d/a),-0.798+ln(d/rw)      +  0.4341n7 J£ 

v  f  m  (rw)r 

Yf  =  j^y-        (33) 

l/Ym[0.682(d/a)f-0.798  +  ln(d/rw)m1  +  1  +  ^M-1)ln-^ 

where  the  subscripts    f    and    m    denote  field  and  model  respectively. 

Fractional  Flow 

Up  to  the  time  of  breakthrough,   the  fractional  flow  from  the 
unswept  region  (oil  cut)  is  one.     Also  up  to  the  time  of  breakthrough, 
a  plot  of  areal  sweep  versus  displaceable  volumes  injected  yields  a 
straight  line  at  a  slope  of  one.     It  can  also  be  seen  that  a  plot  of  areal 
sweep  versus  elapsed  time,   with  injection  at  a  constant  rate,   will 
yield  a  straight  line  prior  to  breakthrough.     After  breakthrough,   the 

areal  sweep  curve  deviates  from  a  straight  line. 

1  3    20 
Previous  publications      '        have  shown  that  the  fractional  flow 

from  the  producing  well  can  be  obtained  from  the  slope  of  the  areal 

sweep  efficiency  curve  versus  displaceable  volumes  injected.     Since 

elapsed  time  and  displaceable  volumes  are  directly  proportional,   the 

fractional  flow  from  the  unswept  region  can  be  obtained  by  the  ratio 

of  the  slope  of  the  areal  sweep  efficiency  curve  at  the  time  of  interest 

to  the  slope  of  the  straight  line  portion.     Tables  of  fractional  flow  are 

given  with  the  basic  data  plots  in  Appendix  B  and  are  plotted  versus 

areal  sweep  efficiency  in  Appendix  C. 


INTERPRETATION  OF  RESULTS 

Tables  of  basic  data  from  the  model  runs  are  presented  in 
Appendix  B.     Appendix  B  also  includes  curves  of  areal  sweep  efficiency 
versus  elapsed  time.     Each  figure  also  includes  tables  of  displaceable 
volumes  injected  and  fractional  flow  data  as  derived  from  these  curves. 
Plots  of  displaceable  volumes  injected,    fractional  flow,   and  conductance 
ratio  versus  areal  sweep  efficiency  for  the  mobility  ratios  studied  are 
provided  in  Appendix  C. 

Results  of  this  study  were  compared  with  those  obtained  by  the 

13 
X-ray  shadowgraph  technique        and  were  found  to  have  a  maximum 

deviation  of  7%  and  an  average  difference  of  less  than  2%.     Comparison 

of  areal  sweep  efficiencies  at  breakthrough  and  at  95%  fractional  flow 

13 
from  the  swept  region  as  obtained  by  Dyes,    et  al       and  as  obtained  in 

this  study  are  listed  in  Table  I  for  the  direct  line  drive  and  in  Table  II 

for  the  staggered  line  drive. 

In  order  to  make  the  results  of  this  study  useful  for  mobility 
ratios  other  than  those  specifically  used  here,    Figures  10  through  15 
indicate  the  effect  of  mobility  ratio  on  performance. 

In  comparing  the  results  obtained  for  the  direct  line  drive    and 
the  staggered  line  drive,    Figure  16  shows  that  areal  sweep  efficiencies 
are  significantly  higher  for  the  staggered  line  drive  at  breakthrough. 
However,    at  a  fractional  flow  from  the  swept  region  of  95  per  cent 
(assumed  economic  limit)  the  areal  sweep  efficiencies  are  similar  for 
the  two  patterns  at  abandonment  up  to  a  mobility  ratio  of  approximately 
10.     At  mobility  ratios  above   10,    previous  studies  have  shown  that  the 

40 
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areal  sweep  efficiencies  decrease  rapidly  in  this  region  for  a  given 
fractional  flow. 
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TABLE  I 

COMPARISON  OF  AREAL  SWEEP  EFFICIENCIES 
FOR  DIRECT  LINE  DRIVE 


Breakthrough 
Sweep  Efficiency 

95$ 

i  Fractional 
Sweep  Effi 

Water  Flow 
.ciency 

Mobility- 
Ratio 

Dy 

ss,    et  a] 

i13 

This  Study 

,13 
Dyes,    et  al 

This  Study 

0.  106 

0.92 

0.85 

0.995 

0.98 

0.274 

0.76 

0.76 

0.99 

0.97 

1.0 

0.55 

0.55 

0.98 

0.96 

3.  13 

0.46 

0.46 

0.97 

0.95 

12.1 

* 

0.41 

# 

0.77 

*Not  available 
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TABLE  II 

COMPARISON  OF  AREAL  SWEEP  EFFICIENCIES 
FOR  STAGGERED  LINE  DRIVE 


Breakthrough 
Sweep  Efficiency 

95% 

Fractional  Water  Flow 
Sweep  Efficiency 

Mobility- 
Ratio 

Dy 

es,    et 

i13 
al 

This  Study 

Dyes,    et  al 

13 

This  Study 

0.106 

0.91 

0.89 

0.995 

0.99 

0.274 

0.85 

0.87 

0.99 

0.98 

1.0 

0.75 

0.77 

0.98 

0.98 

3.13 

0.61 

0.66 

0.95 

0.96 

12.  1 

* 

0.55 

# 

0.89 

*Not  available 
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CONCLUSIONS 

The  results  of  this  study  of  line  drive  flood  patterns  indicate 
that  mobility  ratio  has  a  definite  effect  on  the  areal  sweep  efficiencies 
at  time  of  breakthrough.     At  abandonment  conditions,   the  effect  of 
mobility  ratio  on  areal  sweep  is  observable,    but  not  to  the  pronounced 
extent  as  at  breakthrough  conditions. 

The  conductance  ratios  show  that  the  majority  of  the  flow 
potential  change  takes  place  when  the  front  is  near  the  injection  well 
for  favorable  mobility  ratios  and  near  the  production  well  for  unfavor- 
able mobility  ratios.     At  time  of  breakthrough,    an  apparent  change  in 
conductance  ratio  is  noted. 

The  staggered  line  drive  gave  breakthrough  sweep  efficiencies 
from  4  to  14  per  cent  higher  than  those  obtained  by  the  direct  line  drive. 
At  abandonment  conditions,    only  small  differences  existed  in  areal 
sweep  efficiencies  for  the  direct  and  staggered  line  drive.     However, 
for  very  poor  mobility  ratios,    over  one  additional  displaceable  volume 
of  injected  fluid  was  required  to  reach  abandonment  conditions  with  the 
direct  line  drive. 

The  method  described  for  predicting  performance  of  stratified 
reservoirs  can  be  used  for  any  well  pattern  for  which  conductance 

ratios  and  areal  sweep  efficiencies  are  available.     This  data  is  avail- 

7  20 

able  for  the  five-spot   ,   the  inverted  skewed  seven- spot      ,    and  the 

34 
inverted  nine-spot       patterns.     This  study  allows  similar  predictions 

for  the  direct  and  staggered  line  drive  patterns. 
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The  effect  of  a  free  gas  saturation  on  areal  sweep  efficiency 
and  conductance  ratio  has  not  been  sufficiently  studied.     Accurate 
performance  prediction  for  stratified  reservoirs  requires  additional 
study  on  the  effect  of  free  gas  saturations. 


NOMENCLATURE 

a  well  spacing 

c  constant 

d  row  spacing  of  wells 

Es  fractional  area  swept 

fs  fractional  flow  from  swept  region 

fw  fractional  water  flow 

h  height  or  thickness 

k  permeability 

L  length 

I  partial  length 

M  mobility  ratio 

Np  fractional  oil  produced 

^P  pressure  drop 

q  total  flow 

qw  water  flow 

rf  frontal  radius 

ri  radius  of  radial  flow  region  around  injection  well 

rp  radius  of  radial  flow  region  around  production  well 

rw  radius  of  well  bore 

Sew  saturation  of  connate  water 

Sro  saturation  of  residual  oil 

t  time 

Vb  bulk  volume 

Vd  displaceable  volume  -  Vb  $  (1   -  Sew  -  Sro) 
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v  velocity 

Wi  water  injected 

w  width 

a  angle  open  to  flow  into  wellbore  from  swept  region 

y  conductance  ratio 

$  porosity 

\  mobility  -  k/p. 

g.  viscosity 

Subscripts 

A  radial  flow  area  around  injection  well 

B  linear  flow  area  between  radial  flow  areas 

C  radial  flow  area  around  production  well 

BT  breakthrough 

f  field 

i  initial 

j  any  stratum 

m  model 

o  oil 

s  swept  region 

t  time 

u  unswept  region 

w  water 
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TABLE  III 
DIRECT  LINE  DRIVE 


Run 
4 

Initial  Pressure 
1.7 

Mobility  Ratio 
0.106 

Elapsed 
Time 

Corrected 
Time 

Fractional 
Area 
Swept 

Pressure 

Conductance 
Ratio 

0 

_ 

0.025 

8.6 

0.198 

1.0 

- 

0.057 

9.8 

0.174 

2.0 

- 

0.094 

10.  1 

0.168 

3.0 

- 

0.135 

10.  2 

0.166 

5.0 

- 

0.181 

10.3 

0.165 

19.27 

- 

0.672 

11.2 

0.152 

21.27 

- 

0.741 

11.6 

0.146 

24.27 

- 

0.809 

12.9 

0.132 

27.27 

- 

0.  855 

13.6 

0.125 

29.27 

- 

0.915 

14.5 

0.117 

31.27 

- 

0.937 

14.9 

0.  114 

37.27 

- 

0.954 

15.3 

0.111 

Note:    Areal  sweep  data  is  invalid. 
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TABLE  IV 
DIRECT  LINE  DRIVE 


Run 
9 

Initial  Pressure 
2.2 

Mobility  Ratio 
0.  106 

Elapsed 
Time 

Correcte 
Time 

id         Fractional 
Area 
Swept 

Pressure 

Conductance 
Ratio 

0 

_ 

0.039 

_ 

_ 

1.00 

- 

0.  104 

10.4 

0.216 

2.16 

- 

0.  143 

10.6 

0.208 

9.75 

- 

0.340 

10.8 

0.204 

13.50 

- 

0.479 

10.9 

0.202 

20.08 

- 

0.  712 

11.4 

0.193 

23.41 

- 

0.844 

11.9 

0.185 

24.41 

- 

0.870 

12.6 

0.  175 

25.41 

- 

- 

13.5 

0.163 

26.41 

- 

0.926 

14.2 

0.155 

27.41 

- 

- 

14.7 

0.150 

29.08 

- 

0.949 

15.  3 

0.144 

35.91 

_ 

0.974 

15.4 

0.143 
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TABLE  V 
DIRECT  LINE  DRIVE 


Run 
2 

Init: 

Lai  Pressure 
6.5 

Mobility  Rati 
0.274 

LO 

Elapsed 
Time 

Corrected 
Time 

Fractional 
Area 
Swept 

Pressure 

Conductance 
Ratio 

0 

_ 

0.107 

10.7 

0.607 

3.33 

- 

0.277 

11.1 

0.586 

6.67 

- 

0.453 

11.5 

0.565 

10.00 

- 

0.697 

12.5 

0.520 

13.25 

- 

0.860 

14.1 

0.461 

16.58 

- 

0.925 

14.9 

0.436 

23.08 

- 

0.945 

15.0 

0.433 

26.33 

- 

- 

15.3 

0.425 

29.67 

_ 

0.970 

15.4 

0.422 

Note:    Rerun  as  run  22. 
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TABLE  VI 

DIRECT  LINE  DRIVE 


Run 
9 

Init 

Lai  Pressure 
4.8 

Mobil 
0. 

ity  Rati 
274 

o 

Elapsed 
Time 

Corrected 
Time 

Fractional 
Area 
Swept 

Pressure 

Conductance 
Ratio 

0 

0 

0.194 

7.3 

0.658 

1.0 

1.9 

0.288 

7.3 

0.658 

2.0 

3.8 

0.372 

7.3 

0.658 

3.0 

5.  7 

0.467 

7.3 

0.658 

4.0 

7.6 

0.558 

7.3 

0.658 

5.0 

9.5 

0.665 

7.3 

0.658 

7.0 

13.3 

0.848 

8.5 

0.565 

8.0 

15.2 

- 

9.0 

0.533 

9.0 

17.1 

0.897 

9.2 

0.522 

10.0 

19.0 

- 

9.5 

0.505 

11.0 

20.9 

0.930 

_ 
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TABLE  VII 

DIRECT  LINE  DRIVE 


Run 

1 

Mobility 
1.0 

Rat: 

LO 

Elapsed 
Time 

Corrected 
Time 

Fractional 
Area 
Swept 

Pressure 

C< 

Dnductance 
Ratio 

0 

_ 

0.0944 

1.0 

1.0 

- 

0.240 

C 

1.0 

2.0 

- 

0.385 

O 

1.0 

3.0 

- 

0.519 

N 

1.0 

4.0 

- 

0.639 

S 

1.0 

5.0 

- 

0.  726 

T 

1.0 

6.33 

- 

0.791 

A 

1.0 

7.33 

- 

0.840 

N 

1.0 

10.58 

- 

0.939 

T 

1.0 

17,08 

- 

0.973 

1.0 

Note:    Rerun  as  run  21. 
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TABLE  VIII 

DIRECT  LINE  DRIVE 


Run 
6 

Mobility 
1.0 

Ratio 

Elapsed 
Time 

Corrected 
Time 

Fractional 
Area 
Swept 

Pressure 

Conductance 
Ratio 

0 

0 

0.271 

1.0 

1.0 

0.84 

0.356 

C 

1.0 

2.0 

1.68 

0.450 

O 

1.0 

3.0 

2.52 

0.546 

N 

1.0 

4.0 

3.36 

0.605 

S 

1.0 

5.0 

4.21 

0.667 

T 

1.0 

6.0 

5.05 

0.724 

A 

1.0 

8.0 

6.  72 

0.790 

N 

1.0 

10.0 

8.42 

0.845 

T 

1.0 

15,59 

13.10 

0.934 

1.0 

18.59 

15.60 

0.960 

1.0 
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TABLE  IX 
DIRECT  LINE  DRIVE 


Run 
3 

Initial  Pressure 
13.3 

Mobility 
3.13 

Ratio 

Elapsed 
Time 

Corrected 
Time 

Fractional 
Area 
Swept 

Pressure 

Conductance 
Ratio 

0 

«. 

0.099 

.» 

_ 

1.0 

- 

0.210 

- 

- 

2.0 

- 

0.327 

9.1 

1.46 

3.0 

- 

0.432 

9.1 

1.46 

4.0 

- 

0.482 

8.1 

1.64 

5.0 

- 

0.530 

7.3 

1.83 

10.78 

- 

0.654 

7.0 

1.90 

17.40 

- 

0.750 

7.0 

1.90 

21.25 

- 

0.798 

7.0 

1.90 

28.00 

- 

0.862 

7.0 

1.90 

34.00 

- 

0.896 

7.0 

1.90 

35.00 

_ 

0.913 

7.0 

1.90 

Note:    Areal  sweep  data  is  invalid. 
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TABLE  X 
DIRECT  LINE  DRIVE 


Run 
8 

Initial  Pressure 
12.4 

Mobility 
3.13 

Ratio 

Elapsed 
Time 

Corrected 
Time 

Fractional 
Area 
Swept 

Pressure 

Conductance 
Ratio 

0 

— 

0.051 

_ 

_ 

1.0 

- 

0.179 

9. 

4 

1.32 

2.0 

- 

0.358 

8. 

5 

1.46 

3.0 

- 

- 

8. 

1 

1.53 

4.0 

- 

0.580 

7. 

0 

1.77 

5.0 

- 

0.642 

6. 

5 

1.91 

6.0 

- 

- 

6. 

3 

1.97 

7.0 

- 

0.720 

6. 

2 

2.00 

12.08 

- 

0.829 

6. 

0 

2.07 

18.58 

- 

0.920 

6. 

0 

2.07 

22.42 

- 

0.952 

6. 

0 

2.07 

25.50 

_ 

0.965 

6. 

0 

2.07 
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TABLE  XI 
DIRECT  LINE  DRIVE 


Run 

5 

Initial  Pressure 
14.7 

Mobility 
12.1 

Ratio 

Elapsed 
Time 

Correct* 
Time 

id         Fractional 
Area 
Swept 

Pressure 

Conductance 
Ratio 

0 

— 

0.065 

_ 

_ 

1.0 

- 

0.101 

- 

- 

2.0 

- 

0.  157 

6.7 

2.19 

3.0 

- 

0.207 

6.7 

2.19 

4.0 

- 

0.259 

6.7 

2.19 

6.0 

- 

0.354 

6.7 

2.19 

7.0 

- 

0.410 

6.6 

2.23 

8.0 

- 

- 

5.9 

2.49 

9.0 

- 

0.496 

5.0 

2.94 

10.0 

- 

- 

4.4 

3.34 

11.0 

- 

0.525 

3.5 

4.20 

.12.0 

- 

- 

2.9 

5.07 

13.0 

- 

0.545 

2.8 

5.26 

18.62 

- 

0.586 

2.7 

5.45 

22.93 

- 

- 

2.6 

5.65 

29.13 

- 

0.616 

2.3 

6.40 

39.67 

- 

0.670 

2.3 

6.40 

56.93 

- 

0.709 

2.3 

6.40 

74.25 

- 

0.742 

2.3 

6.40 

91.55 

- 

0.778 

2.3 

6.40 

119.38 

_ 

0.812 

2.3 

6.40 

Note:    Areal  sweep  data  is  invalid. 
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TABLE  XII 
DIRECT  LINE  DRIVE 


Run 
10 

Initial  Pressure 
15.0 

Mobility 
12.1 

Ratio 

Elapsed 
Time 

Correcte 
Time 

id         Fractional 
Area 
Swept 

Pressure 

Conductance 
Ratio 

0 

_ 

0.038 

_ 

_ 

1.0 

- 

0.107 

9.1 

1.65 

2.0 

- 

0.173 

9.1 

1.65 

3.0 

- 

0.247 

8.8 

1.70 

4.0 

- 

0.323 

8.5 

1.77 

5.0 

- 

0.401 

8.3 

1.81 

6.0 

- 

- 

7.0 

2.14 

7.0 

- 

0.507 

4.2 

3.57 

8.0 

- 

- 

3.6 

4.17 

9.0 

- 

0.551 

3.4 

4.41 

10.0 

- 

- 

3.4 

4.41 

18.0 

- 

0.654 

3.0 

5.00 

24.59 

- 

0.706 

2.8 

5.36 

28.34 

- 

- 

2.8 

5.36 

34.92 

- 

0.740 

2.8 

5.36 

41.50 

- 

0.789 

2.8 

5.36 

51.92 

- 

0.807 

2.8 

5.36 

62.25 

_ 

0.833 

2.8 

5.36 
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TABLE  XIII 

STAGGERED  LINE  DRIVE 


Run 
14 

Initial  Pressure 
2.7 

Mobility  Ratio 
0.106 

Elapsed 
Time 

Corrected         Fractional 
Time                   Area 
Swept 

Pressure 

Conductance 
Ratio 

0 

_ 

0.122 

11.7 

0.230 

6.08 

- 

0.282 

11.8 

0.229 

13.08 

- 

0.482 

12.0 

0.225 

16.66 

- 

0.598 

12.1 

0.223 

23.16 

- 

0.  768 

12.5 

0.216 

26.16 

- 

0.918 

13.7 

0.197 

27.16 

- 

- 

14.2 

0.190 

28.16 

- 

0.955 

14.7 

0.184 

29.16 

- 

0.969 

15.0 

0.180 

30.16 

- 

- 

15.2 

0.178 

31.16 

- 

- 

15.4 

0.175 

32.16 

- 

0.978 

15.6 

0.173 

37.16 

_ 

0.989 

15.7 

0.171 

TABLE  XIV 
STAGGERED  LINE  DRIVE 


Run 

19 

Initial  Pressure 
2.8 

Mobility  Ratio 
0.  106 

Elapsed 
Time 

Corrected 
Time 

Fractional 
Area 
Swept 

Pressure 

Conductance 
Ratio 

0 

0 

0.0086 

9.8 

0.284 

1.87 

0.505 

0.032 

12.1 

0.231 

5.20 

2.21 

0.074 

12.5 

0.224 

12.12 

5.17 

0.  174 

12.  5 

0.224 

18.62 

7.94 

0.279 

12.5 

0.224 

25.37 

10.80 

0.384 

12.5 

0.224 

32.03 

13.65 

0.503 

12.5 

0.224 

38.62 

16.45 

0.630 

12.  5 

0.224 

45.20 

19.20 

0.717 

12.5 

0.224 

51.78 

22.  10 

0.822 

12.5 

0.224 

55.12 

23.50 

0.867 

13.0 

0.215 

58.45 

24.90 

0.900 

14.0 

0.200 

61.70 

26.30 

0.932 

15.0 

0.187 

64.95 

27.60 

0.948 

15.8 

0.177 

74.70 

31.80 

0.975 

16.3 

0.172 

84.37 

35.90 

0.990 

16.8 

0.170 

82 


1 

o  r^  r\j  so  a^  ^  ^h 
O  00  (M  co  co  ^  00 
o  -*  iTi  r*-  oo  o  o 

o  o  o  o  o  o  o 

T3 
> 

o  o  o  o  o  o  o 
o  <\j  o  vO  co  o  r* 

00  O  O  O  —i  f\]  CM 

o  O  O  •—*  •—*  ^h  i— < 

- 

w 

OOiTOO^N 
O  >-<   ^   Is-  00  00  O 
00  O^  O^  0^  0s  O^  0s 

o  o  o  o  o  o  o 

•H 

H 

-<  o  o  ©  o  o  © 

in  o  oo  o  cj  i*  o 

CM  CM  c\]  CO  CO  CO  CO 

- 

K^P 

- 

<* 

l—t 

>v     X 

_ 

a 

03 

X 

0 

1             1 

1                 1                 1                  1                 1                  1 

c>  K 

o 


cm 

rO 


_     "•£> 


-    oo 


00 

o 


o 


Efl 


O 


O 


W 
2 


Q 
W 

CO 
Ph 

W 

CO 

CO 

U 
> 


CO 

03 

< 

i 

i 

w 
> 

I— I 

03 

Q 


g  w 

6  s 

S'3 

Q 
W 

03 

W 

a 
a 
< 

H 
co 


r- 
N 

w 
o 

0 

M 


83 


TABLE  XV 
STAGGERED  LINE  DRIVE 


Run 

Initial  Pressure 

Mobility  Ratio 

12 

5.5 

0.274 

Elapsed 

Corrected         Fractional 

Pressure 

Conductance 

Time 

Time                   Area 
Swept 

Ratio 

0 

0.066 

- 

- 

2.0 

0.187 

- 

- 

4.0 

0.275 

— 

- 

5.0 

0.347 

- 

- 

6.0 

0.396 

10.7 

0.514 

7.0 

0.461 

10.7 

0.514 

9.0 

0.584 

10.7 

0.514 

11.0 

0.718 

10.7 

0.514 

13.0 

0.842 

11.3 

0.486 

14.0 

0.923 

11.8 

0.466 

15.0 

0.950 

12.0 

0.458 

17.0 

0.971 

12.2 

0.451 

21.0 

0.980 

12.2 

0.451 

25.0 

0.988 

12.2 

0.451 
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TABLE  XVI 

STAGGERED  LINE  DRIVE 


Run 
17 

Initi 

Lai  Pressure 
5.1 

Mobility  Ratio 
0.274 

Elapsed 
Time 

Corrected 
Time 

Fractional 
Area 
Swept 

Pressure 

Conductance 
Ratio 

0 

0 

0.414 

9.1 

0.560 

1.0 

0.934 

0.871 

9.5 

0.537 

2.0 

1.88 

0.143 

9.5 

0.537 

3.0 

2.80 

0.196 

9.6 

0.531 

4.0 

3.74 

0.257 

9.6 

0.531 

6.0 

5.60 

0.372 

9.6 

0.531 

8.0 

7.47 

0.481 

9.6 

0.531 

10.0 

9.34 

0.597 

9.6 

0.531 

12.0 

11.20 

0.703 

9.7 

0.526 

13.0 

12.10 

- 

9.8 

0.520 

14.0 

13.10 

0.793 

10.1 

0.505 

15.0 

14.00 

0.884 

10.6 

0.481 

16.0 
17.0 

14.94 
15.88 

0.951 

11.0 
11.  3 

0.464 
0.451 

19.0 

17.75 

0.976 

11.5 

0.443 

20.0 

18.68 

- 

11.6 

0.440 

21.0 

19.60 

0.987 

11.7 

0.436 
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TABLE  XVII 

STAGGERED  LINE  DRIVE 


Run 
11 

Mobility 
1.0 

Ratio 

i 

Elapsed 
Time 

Corrected 
Time 

Fractional 
Area 
Swept 

Pressure 

Conductance 
Ratio 

0 

0 

0.083 

C 

1.0 

1.0 

1.05 

0.272 

O 

1.0 

2.0 

2.10 

0.473 

N 

1.0 

3.0 

3.15 

0.685 

S 

1.0 

4.0 

4.20 

0.823 

T 

1.0 

5.0 

5.25 

0.904 

A 

1.0 

6.0 

6.30 

0.951 

N 

1.0 

8.0 

8.45 

0.987 

T 

1.0 
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TABLE  XVHI 
STAGGERED  LINE  DRIVE 


Run 
16 

Mol 

>ility 
1.0 

Ratio 

Elapsed 
Time 

Corrected 
Time 

Fractional 
Area 
Swept 

Pressure 

Conductance 
Ratio 

0 

. 

0.113 

C 

1.0 

1.0 

- 

0.281 

O 

1.0 

2.0 

- 

0.450 

N 

1.0 

3.0 

- 

0.611 

S 

1.0 

4.0 

- 

0.775 

T 

1.0 

5.0 

- 

0.876 

A 

1.0 

6.0 

- 

0.940 

N 

1.0 

7.0 

— 

0.967 

S 

1.0 

8.0 

»— 

0.992 

T 

1.0 
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TABLE  XIX 

STAGGERED  LINE  DRIVE 


Run 
13 

Initial  Pressure 
14.7 

Mobility 
3.13 

Ratio 

Elapsed 
Time 

Corrected         Fractional 
Time                   Area 
Swept 

Pressure 

Conductance 
Ratio 

0 

_ 

0.098 

9.1 

1.62 

1.0 

- 

0.229 

8.9 

1.65 

2.0 

- 

0.369 

8.9 

1.65 

3.0 

- 

0.509 

8.8 

1.67 

4.0 

- 

0.636 

8.8 

1.67 

5.0 

" 

0.739 

8.1 

1.81 

6.0 

— 

0.778 

7.3 

2.01 

7.0 

— 

- 

7.1 

2.07 

8.0 

- 

0.832 

7.0 

2.10 

10.0 

- 

0.873 

7.0 

2.10 

13.0 

- 

0.911 

7.0 

2.10 

15.0 

- 

0.936 

7.0 

2.10 

17.0 

- 

0.969 

7.0 

2.10 

19.0 

- 

0.988 

7.0 

2.10 
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TABLE  XX 

STAGGERED  LINE  DRIVE 


Run 
18 

Initial  Pressure 
14.2 

Mobility 
3.  13 

Ratio 

Elapsed 
Time 

Corrected         Fractional 
Time                   Area 
Swept 

Pressure 

Conductance 
Ratio 

0 

0 

0.098 

_ 

_ 

1.0 

0.728 

0.191 

- 

- 

2.0 

1.46 

0.287 

- 

- 

3.0 

2.19 

0.377 

- 

- 

4.0 

2.92 

0.465 

- 

- 

5.0 

3.64 

0.553 

10.5 

1.35 

6.0 

4.37 

0.632 

10.5 

1.35 

7.0 

5.  10 

0.694 

10.4 

1.37 

8.0 

5.83 

0.759 

10.  1 

1.41 

9.0 

6.56 

0.802 

9.9 

1.43 

10.0 

7.28 

- 

9.5 

1.49 

11.0 

8.02 

0.837 

9.2 

1.54 

12.0 

8.75 

- 

9.1 

1.56 

13.0 

9.48 

0.866 

9.0 

1.58 

15.0 

10.92 

0.895 

9.0 

1.58 

18.0 

13.10 

0.920 

9.0 

1.58 

20.0 

14.60 

0.956 

9.0 

1.58 

23.0 

16.80 

1.973 

9.0 

1.58 
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TABLE  XXI 

STAGGERED  LINE  DRIVE 


Run 

15 

Init 

Lai  Pressure 
16.0 

Mobility 
12.1 

Ratio 

Elapsed 
Time 

Corrected 
Time 

Fractional 
Area 
Swept 

Pressure 

Conductance 
Ratio 

0 

_ 

0.014 

_ 

- 

1.0 

- 

0.072 

8.6 

1.86 

2.0 

- 

0.134 

8.2 

1.95 

3.0 

- 

- 

7.8 

2.05 

4.0 

- 

0.254 

7.6 

2.11 

5.0 

- 

0.315 

7.4 

2.16 

7.0 

— 

0.424 

7.1 

2.25 

8.0 

- 

0.477 

7.0 

2.29 

9.0 

- 

- 

6.3 

- 

2.54 

10.0 

- 

0.584 

4.8 

3.33 

12.0 

- 

0.644 

4.8 

3.33 

13.0 

- 

- 

4.2 

3.81 

14.0 

- 

- 

3.8 

4.21 

16.0 

- 

0.  715 

3.4 

4.71 

18.0 

- 

0.749 

3.4 

4.71 

24.91 

- 

0.785 

3.4 

4.71 

42.65 

- 

0.877 

3.4 

4.71 

55.08 

_ 

0.915 

3.4 

4.71 
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TABLE  XXII 

STAGGERED  LINE  DRIVE 


Run 

Initial  Pressure 

Mobility 

Ratio 

20 

16.5 

12.1 

Elapsed 

Corrected 

Fractional 

Pressure 

Conductance 

Time 

Time 

Area 
Swept 

Ratio 

0 

0 

0.038 

9.2 

1.79 

1.0 

0.523 

0.  104 

8.2 

2.01 

4.33 

2.27 

0.204 

7.8 

2.11 

7.58 

3.97 

0.294 

7.8 

2.11 

10.91 

6.79 

0.385 

7.8 

2.11 

14.16 

7.42 

0.468 

7.8 

2.11 

17.50 

9.15 

0.553 

7.8 

2.11 

18.50 

9.68 

0.579 

7.8 

2.11 

20.50 

10.  7 

- 

7.5 

2.20 

21.50 

11.2 

0.648 

7.2 

2.29 

22.50 

11.8 

- 

7.0 

2.36 

25.00 

13.1 

0.693 

6.8 

2.43 

28.25 

14.8 

- 

6.2 

2.66 

31.66 

16.6 

0.747 

5.5 

3.0 

39.25 

20.5 

0.771 

5.5 

3.0 

48.08 

25.  1 

0.807 

5.5 

3.0 

57.83 

30.2 

0.827 

5.5 

3.0 

67.41 

35.2 

0.853 

5.5 

3.0 

76.08 

39.8 

0.872 

5.5 

3.0 
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APPENDIX  C 


FIGURES  OF  INTERPRETED  DATA 
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FIGURE  32.     DIRECT  LINE  DRIVE  --  MOBILITY  RATIOS 

for  DISPLACEABLE  VOLUMES  INJECTED  versus 
AREA  SWEPT 
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FIGURE  33.     DIRECT  LINE  DRIVE  --  MOBILITY  RATIOS 

for  FRACTIONAL  FLOW  versus  AREA  SWEPT 
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FIGURE  34.     DIRECT  LINE  DRIVE  --  MOBILITY  RATIOS 

for  CONDUCTANCE  RATIO  versus  AREA  SWEPT 
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FIGURE  35.  STAGGERED  LINE  DRIVE  --  MOBILITY  RATIOS 
for  DISPLA.CEABLE  VOLUME  INJECTED  versus 
AREA  SWEPT 
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FIGURE  36.     STAGGERED  LINE  DRIVE  --  MOBILITY  RATIOS 
for  FRACTIONAL  FLOW  versus  AREA  SWEPT 
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FIGURE  37.     STAGGERED  LINE  DRIVE  --  MOBILITY  RATIOS 
for  CONDUCTANCE  RATIO  versus  AREA  SWEPT 


APPENDIX  D 


COMPUTED  CONDUCTANCE  RATIOS 
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PROGRAM  CONRA 

DIMENSION  RM( 10) ♦  DFS (20) »ES( 10*20) 

READ  1.  JE»(RMU)»  J=1»JE) 

READ  It  KE,  (DFS(K),  K=1,KE) 

1  FORMAT  (  I3,11F7.0) 
DO  10  J=1,JE 

2  FORMAT  (11F7.0) 

10  READ  2»  (ES(J>K)»  K=1,KE) 

RM  IS  THE  MOBILITY  RATIO  AND  JE  IS  THE  NUMBER  OF  MOBILITY  RATIOS 
WHICH  ARE  TO  BE  STUDIED.   DFS  IS  FRACTIONAL  FLOW  FROM  THE  SWEPT  REGION 
AND  KE  IS  THE  NUMBER  OF  FRACTIONAL  FLOW  VALUES  TO  BE  READ  IN. 
ES  VALUES  ARE  AREAL  SWEEP  EFFICIENCIES  CORRESPONDING  TO  THE  DFS  VALUES. 
SET  INITIAL  VALUES 
D  =  1.0 
A  =  D 

REI  =  0.25*D 
RW  =  D/320.0 

RLW=  0.5*(D/A  -  1.17  +  0.637*LOGF ( A/REI ) ) 
ESA  =  3.1416*REI**2/(2.0*D*A) 
DO  70  J=1,JE 
PRINT  3,RM(J) 
30FORMAT  ( 1H1 ,40X , 16HMOB I L I TY  RATIO  =F5 . 2 // 10X , 19H  CONDUCTIVITY  RATI 
10, 10X»  INFRACTION  SWEPT  *  1  3X  ♦  2HFS//2  3X  ,  7H1  .  00000  » 19  X  ,  1H0  »  19X  » 1H0  ) 
ESBT  =  ES(J,2) 
REP  =  SQRTF(2.0*D*A*( 1 . O-ESBT ) /3 . 1416 ) 


104 

DES  =  0.01 

ES(J»1)  =  DES 

DO  70  K»1»KE 
80  IF  (ES(J»K)  -  ESA)  20,20*30 
30  IF  (ES(J,K)  -  ESBT)  40,50,50 
20  DLW  =  0.0 

RF  ■  SQRTF(2.0*D*A*ES( J,K)/3.1416) 

GO  TO  60 
40  RF  =  REI 

DLW  =RLW*(ES( J,K)  -  ESA)/(ESBT  -  ESA) 

DES  =  0.05 

GO  TO  60 
50  RF  =  REI 

DLW  =  RLW 
600CR  =  (LOGF(REI/RW)  +  6.28*RLW+  LOGF ( REP/RW ) ) / ( 1 . /RM( J )*LOGF ( RF/RW ) 
1+  LOGF(REI/RF)  +  6.28* ( 1 . /RM( J )  -  1.0)*DLW  +  6.28*RLW+  LOGF(REP/RW 
2)*(1.0  -  DFS(K) ) /(1.0-DFS(K)/(RM( J)-DFS(K)*(RM( J)-1.0) ) ) ) 

PRINT  4.  CR,ES( J»K),DFS(K) 
4  FORMAT  (10X.3F20.5) 

ES(J,K)  =  ES(J,K)  +  DES 

IF  (ES(J,K)  -  ESBT)  80,70,70 
70  CONTINUE 

END 
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Flow  Diagram  Program  CONRA 


Dimension 
M,  fs,   Es 


Read 

M,   fs,   Es 


Set  initial  values 

d  =  1.0 

a  =  d 

ri  =  0.25d 

rw  =  d/  320.0 

L/w  =  0.5(d/a  -  1.  17  +  0.  637  loga/ri  ) 

Es.    =  n  ri^  /  2da 


Do  Loop 
for  all 
M  values 


EsBT  =  Es(j,2) 

rp  =  (2da(l  -EsBT)  /  n) 

AEs  =  0.01 

Es(j,l)  =  AEs 


0.5 


Do  Loop 
for  all 
fs  values 
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Yes 


i/v/  =  0 

fs  =  0 

rf  =  (2daEs  /tt) 


0.  5 


No 


Yes 


No 


rf  =  ri 

fs  =  0 

AEs  =  0.05 


X/w 


_L(Es  -  EsA) 
w(EsBT'  ESA] 


rf  =  ri 
i/w  =  L/w 


1 


Y  =  (lnri/rw  +  6.28L/w  +  lnrp/rw)  / 

<((l/M)lnrf/rw  +  lnri/rf  +  6.28(1/M  -   1)  4/w 
+  6.28L/w  +  lnrp/rw    |j  -  f  s  /  (M  -  fs(M  -  1))]^> 


Print 

Y,   Es,   fs 


Yes 


No 
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FORTRAN 

Nomenclature 

A 

a 

CR 

Y 

D 

d 

DES 

AEs 

ES 

Es 

ESA 

ESA 

ESBT 

ESBT 

DFS 

fs 

DLW 

X/w 

REI 

ri 

Description 


REP 


rp 


well  spacing 

conductance  ratio 

row  spacing  of  wells 

increment  of  area  swept 

fractional  area  swept 

radial  flow  region  around 
injection  well 

area  swept  at  breakthrough 

fractional  flow  from  swept  region 

partial  length  to  width  ratio 

radius  of  radial  flow  region 
around  injection  well 

radius  of  radial  flow  region 
around  production  well 


RF 

rf 

frontal  radius 

RLW 

L/w 

length  to  width  ratio 

RM 

M 

mobility  ratio 

RW 

rw 

radius  of  wellbore 
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FIGURE  38.     DIRECT  LINE  DRIVE  --  COMPARISON  OF 

ANALYTICAL  and  MODEL  CONDUCTANCE  RATIOS 
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FIGURE  39.     STAGGERED  LINE  DRIVE  --  COMPARISON  OF 

ANALYTICAL  and  MODEL  CONDUCTANCE  RATIOS 
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FIGURE  40.     DIRECT  LINE  DRIVE  --  EFFECT  OF  L/w  RATIO 
ON  CONDUCTANCE  RATIO 
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APPENDIX  E 
EXAMPLE  PREDICTION  FOR  STRATIFIED  RESERVOIR 

To  illustrate  the  procedure  described  for  stratified  reservoirs, 
assume  a  staggered  line  drive  unit  of  40  acres  well  spacing  composed 
of  2  strata.     Under  reservoir  conditions,   the  oil  viscosity  is  4  cp  and 
the  water  viscosity  is  0.  4  cp.     Assume  that  no  gas  saturation  exists. 

The  reservoir  properties  are: 

Height 

Porosity 

Connate  water  (Sew) 

Residual  oil  (Sro) 

k     (at  Sew) 

o  x  ' 

k     (at  Sro) 

w  x  ' 

The  bulk  volume  in  each  stratum  will  be  the  same  or  (20)  (80)  = 

1600  acre  ft. 

Vd      =    Vb  $  (1-Scw-Sro) 

Vdx    =    (1600)  (0.30)  (1-0.30-0.30)    =    192  acre-ft 

Vd2    =    (1600)  (0.  25)  (1-0.25-0.30)    =    180  acre-ft 

The  mobility  ratios  for  each  stratum  are 

m   =  (Viv)/(Vn.) 

Mj     =    (20/0.4)/(200/4)  =  1 


Stratum  1 

Stratum  2 

20  feet 

20  feet 

0.30 

0.25 

0.30 

0.25 

0.30 

0.30 

200  md 

75  md 

20  md 

22  md 

M 


2    =    (22/0.4)/(75/4)  =  3.13 


and 
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113 
ETW;    =    175)  (180)    =    z'84 


klVdl  (200) (192) 


^2Y~2 

To  obtain    ct    for  each  stratum  from  equation  (25)  the  inverse 
of  the  conductance  ratio  is  plotted  versus  displaceable  volumes  injected 
(from  Figures  13  and  15).     The  resulting  curves  (Figure  41)  are 
graphically  integrated  to  obtain    ct,     and    ct-     .     (Note  that  the  data  for 
each  Figure  is  in  the  preceding  Table. ) 

Cumulative  water  injected  and  the  corresponding  oil  produced 
are  obtained  from  Figure  13  in  dimensionless  terms  of  displaceable 
volumes  injected  and  fractional  area  swept.     These  are  then  converted 
to    Wi    and    Np    and  plotted  for  each  stratum  versus    ct    in  Figure  42. 
Total  water  injected  and  oil  produced  can  be  obtained  from  Figure  42 
by  summing  the     Wi    and    Np    curves  for  each  stratum  and  is  plotted 
in  Figure  43.     This  curve  can  be  converted  to  an  actual  time  relation- 
ship by  knowing  the  maximum  injection  rate  and  well  allowables. 
Fractional  flow  from  the  swept  region  is  obtained  as  the  slope  of  the 
Np    versus    Wi    curve  and  the  total  conductance  ratio  is  obtained  from 
equation  (31).     Fractional  flow  is  plotted  on  Figure  43  and  conductance 
ratio  versus  water  injected  for  the  system  is  plotted  on  Figure  44. 

The  above  example  does  not  consider  the  presence  of  a  free  gas 
saturation  at  the  start  of  the  water  flood.     However,   it  will  be  within 
the  realm  of  engineering  accuracy  to  assume  that  oil  bank  fill  up  occurs 
simultaneously  in  all  strata.     This  will  be  true  except  in  the  case  of 
high  gas  saturations  because  of  the  unusually  favorable  mobility  ratio 
during  oil  bank  fill  up.     The  volume  occupied  by  free  gas  prior  to  oil 
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bank  fill  up  can  then  be  added  to  the  cumulative  water  injected  and 

volume  swept  to  obtain  satisfactory  performance  predictions  when  there 

is  free  gas  present  at  the  beginning  of  the  flood. 
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TABLE  XXIII 


Wi/Vd 
fraction 

Wi 
acre -ft 

Y2 

i/y2 

0 

0 

1.0 

1.0 

0.1 

18.0 

1.48 

0.675 

0.5 

90.0 

1.50 

0.675 

0.665 

119.8 

1.52 

0.657 

1.0 

180.0 

1.75 

0.571 

1.5 

270.0 

1.83 

0.546 

2.0 

360.0 

1.84 

0.544 

2.5 

450.0 

1,85 

0.540 
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117 


Stratum  one 

Stratum  two 

Wi 

ct. 

ct2 

Wi 

(1/y)  dWi 

2.84   (1/Y)  dWi 

0.0 

0. 

0 

0.0 

0.0 

0.0 

100.0 

100, 

0 

50.0 

35.  38 

100.  5 

200.0 

200 

0 

100.0 

68.83 

195.5 

300.0 

300. 

0 

150.0 

100.88 

285.0 

400.0 

400. 

0 

200.0 
300.0 
400.0 
500.0 

132.88 
188. 18 
242.48 
296. 38 

377.0 
535.0 
688.0 
843.0 

Wi/Vd       Np/Vd 

Wi 

Np 

Wi/Vd 

Np/Vd 

Wi             Np 

fraction    fraction 

acre -ft 

acre- 

•ft 

fraction 

fraction 

acre -ft    acre- ft 

0.769         0.769 

147.5 

147. 

5 

0.665 

0.665 

119;5        119.5 

0.80           0.795 

153.  5 

152. 

5 

0.80 

0.743 

143.8        133.5 

1.0             0.90 

192.0 

173. 

0 

1.0 

0.803 

180.0        144.2 

1.2             0.952 

230.0 

182. 

8 

1.5 

0.883 

269.0        158.7 

1.4             0.985 

269.0 

189. 

0 

2,0 

2.5 

0.930 
0.961 

359,0        167.0 
449.0        173.0 
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TABLE  XXV 


ct 

Wij 

wi2 

Npj 

Np2 

Wi 

Np 

100 

100 

50 

100 

50 

150 

150 

145 

145 

75 

145 

75 

220 

220 

200 

200 

105 

173 

105 

305 

278 

225 

225 

120 

180 

120 

345 

300 

250 

250 

132 

185 

127 

382 

312 

300 

300 

160 

190 

138 

460 

328 

350 

350 

187 

190 

146 

537 

336 

400 

400 

216 

190 

152 

616 

342 

500 

500 

276 

190 

160 

776 

350 

600 

600 

344 

180 

166 

944 

356 

700 

700 

410 

190 

171 

1110 

361 

Wi 

fo 

fw 

220 

1.000 

0.000 

300 

0.590 

0.410 

340 

0.550 

0.450 

350 

0.400 

0.600 

400 

0.245 

0.755 

500 

0,127 

0.873 

600 

0,054 

0.946 

800 

0.038 

0.962 

1000 

0.031 

0.969 
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TABLE  XXVI 


k2Vd2 


Yl  +  kjVdj    Y2  1  +  0.352  y2 

1  + 


V  k2Vd2 1.352 


r^r 


ct 

wi2 

1/Y2 

Y2 

Wi 

Y 

50 

25 

0.675 

1.480 

75 

1.126 

100 

50 

0.672 

1.490 

150 

1.128 

150 

78 

0.670 

1.492 

225 

1.128 

200 

105 

0.660 

1.518 

305 

1.134 

300 

160 

0.600 

1.669 

460 

1.172 

400 

216 

0.  552 

1.  812 

616 

1.210 

500 

276 

0.550 

1.820 

776 

1.214 

600 

344 

0.542 

1.845 

944 

1.220 

700 

410 

0.540 

1.850 

1110 

1.221 
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